Oscillating Bianchi IX Universe in Hofava-Lifshitz Gravity 
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We study a vacuum Bianchi IX universe in the context of Hoi^ava-Lifshitz (HL) gravity. In 
particular, we focus on the classical dynamics of the universe and analyze how anisotropy changes 
the history of the universe. For small anisotropy, we find an oscillating universe as well as a 
bounce universe just as the case of the Friedmann-Lemaitre-Robertson- Walker (FLRW) spacetime. 
However, if the initial anisotropy is large, we find the universe which ends up with a big crunch 
after oscillations if a cosmological constant A is zero or negative. For A > 0, we find a variety of 
histories of the universe, that is de Sitter expanding universe after oscillations in addition to the 
oscillating solution and the previous big crunch solution. This fate of the universe shows sensitive 
dependence of initial conditions, which is one of the typical properties of a chaotic system. If the 
initial anisotropy is near the upper bound, we find the universe starting from a big bang and ending 
up with a big crunch for A < 0, while de Sitter expanding universe starting from a big bang for 
A > 0. 

PACS numbers: 04.60.-m, 98.80.Cq, 98.80.-k 
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I. INTRODUCTION 



Since the advent of the big bang theory, the initial sin- 
gularity problem is of prime importance in the field of 
cosmology. As shown by Hawking and Penrose [l|, gen- 
eral relativity (GR) predicts a spacetime singularity if a 
certain condition is satisfied. Their singularity theorem 
concludes that our universe must have an initial singu- 
larity. However, once a singularity is formed, general rel- 
ativity becomes no longer valid. It must be replaced by 
more fundamental gravity theory. Even in the framework 
of an inflationary scenario which resolves many difficul- 
ties in the early universe based on the big bang theory, 
the initial singularity cannot be avoided. New gravita- 
tional theory may be required to describe the beginning 
of the universe. 

Many researchers attempt to resolve this singularity 
problem in the context of generalization or extension 
of general relativity^. However, no success has been 
achieved yet. Superstring theory, which is one of the 
most promising candidates for unified theory of funda- 
mental interactions, may solve it, but so far it has not 
been completed yet and is not so far able to describe any 
realistic strong gravitational phenomena. Loop quantum 
gravity theory may resolve the problem of a big bang 
singularity via loop quantum cosmology Q. However it 
is still unclear how to describe time evolution of quantum 
spacetime in loop quantum gravity because of the lack of 
"time" variable. 
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Among attempts to construct a complete quantum 
gravitational theory, Hofava-Lifshitz (HL) gravity has 
been attracted much interest as a candidate for such 
a theory over the past years. HL gravity is charac- 
terized by its power-counting renormalizablity, which is 
brought about by a Lifshitz-like anisotropic scaling as 
t — >■ £^t,x — >■ £x, with the dynamical critical exponent 
z = 3 in the ultra-violet (UV) limit In order to 
recover general relativity (or the Lorentz invariance) in 
our world, one expects that the constant A converges to 
unity in the infrared (IR) limit in the renormalization 
flow. Although it has been argued that there exist some 
fundamental problems in HL gravity [sl-flGj. some exten- 
sions are proposed to remedy these difficulties [Tol. [T7l-[l9j . 
It is intriguing issue whether or not HL gravity can be a 
complete theory of quantum gravity. 

There are a number of works on cosmology in HL grav- 
ity 0, [20l - l53j . As pointed out by earlier works, a big 
bang initial singularity may be avoided in the frame- 
work of HL cosmology due to the higher order terms 
in the spatial curvature %ij in the action [ssj . In this 
context, many researchers have studied the dynamics 
of the Friedmann-Lemaitre-Robertson- Walker (FLRW) 
universe in HL gravity [35l - l5l| . In isotropic and ho- 
mogeneous spacetime, higher curvature terms with ar- 
bitrary coupling constants mimic various types of matter 
with arbitrary sign of energy densities. The 2 = 2 and 
2 = 3 scaling terms give "dark radiation" and "dark stiff- 
matter" , respectively. Although "dark radiation" terms 
in the models with the detailed balance condition can 
avoid the initial singularity, such terms may become irrel- 
evant to the dynamics when we include rclativistic mat- 
ter fields, which may scale as z = 3 in the UV limit and 
behave Stiff matter [U. 

In our previous paper [47| , we have studied the dynam- 
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ics of vacuum FLRW spacetimc in generalized HL gravity 
model without the detailed balance condition and shown 
that "dark stiff-matter" can avoid the initial singular- 
ity of the universe. Even if we include relativistic mat- 
ter fields, when the contribution of "dark stiff-matter" 
is dominant, the singularity is avoided and an oscillating 
spacetime or a bounce universe is obtained. 

Although we have shown a singularity avoidance in HL 
cosmology, the following question may arise: Is this sin- 
gularity avoidance generic? Is such a non-singular space- 
time stable against anisotropic and/or inhomogeneous 
perturbations? In order to answer for these questions, 
we have to study more generic spacetime than the FLRW 
universe. 

The initial state of the universe could be anisotropic 
and/or inhomogeneous. Before the singularity theorem, 
some people believed that the big bang singularity ap- 
pears because of its high symmetry and it may be re- 
solved if one studies anisotropic and/or inhomogeneous 
spacetime. Then they analyzed anisotropic Bianchi- 
type universes and their generalization. Although they 
found some interesting behaviours near the singularity 
such as chaos in Bianchi IX spacetime [54l - l57l |. they 
could not succeed the singularity avoidance. It is sim- 
ply because the singularity theorem does not allow a 
singularity avoidance in GR. The situation becomes 
worse if we consider anisotropy and/or inhomogeneity. 
Even in the effective gravity model derived from super- 
string, which shows a singularity avoidance, with such a 
property [5^ [5^, once we include anisotropy and/or inho- 
mogeneity, the property of such a singularity avoidance 
may be spoiled |6C|. 

Therefore it is important to study whether or not non- 
singular universes in the present HL cosmology still exist 
with anisotropy and/or inhomogeneity. In the present 
paper, we shall investigate the possibility of the singular- 
ity avoidance in homogeneous but anisotropic Bianchi IX 
universe. Since we are interested in a singularity avoid- 
ance, we focus on an oscillating universe and analyze how 
anisotropy changes the history of the universe. We will 
not study the chaotic behaviour in detail, which may ap- 
pear near the big bang singularity, although it is one of 
the most popular and important properties in the Bianchi 
IX spacetime and was discussed analytically in [52], HI] . 
As we will show later, however, some property of non- 
integrable system, i.e., sensitive dependence on initial 
conditions may be found in the fate of the universe in 
the present analysis as well. 

The paper is organized as follows. After a short 
overview of HL gravity, we present the basic equations 
for the vacuum Bianchi IX universe in HL gravity in 
Secini In Sec lIIH we study the stability of the closed 
FLRW universe against small anisotropic perturbations. 
In Sec lIVI we analyze Bianchi IX universe numerically 
and show a variety of histories of the universe, depend- 
ing on initial anisotropy. Summary and remarks follow 
in Sec|Vl In Appendix, we also analyze a bounce uni- 
verse with anisotropy as an another type of non-singular 



solution. 



II. BIANCHI IX UNIVERSE IN 
HORAVA-LIFSHITZ GRAVITY 

First we introduce our Lagrangian of HL gravity, by 
which we will discuss the Bianchi IX universe. The basic 
variables in HL gravity are the lapse function, N, the 
shift vector, Ni, and the spatial metric, gij. These vari- 
ables are subject to the action d, [l^ 

^HL = ^ y" dt£x^N {^K - rHL[ff,,]) , (2.1) 

where = l/Mp^ (Mpl: the Planck mass) and the 
kinetic term is given by 

^K='K^J'K'' -X'K' (2.2) 

with 

■■= ^ (.dtg^j - V^N, - V,iV,) (2.3) 

■K■■^g''yC^J (2.4) 

being the extrinsic curvature and its trace. The poten- 
tial term will be defined shortly. In GR we have 
A = 1, only for which the kinetic term is invariant un- 
der general coordinate transformations. In HL gravity, 
however, Lorentz symmetry is broken in exchange for 
renormalizability and the theory is invariant under the 
foliation-preserving diffeomorphism transformations, 

t^t{t), x'^x'{t,x^). (2.5) 

As implied by the symmetry (|2.5p . it is most natu- 
ral to consider the projectable version of HL gravity, for 
which the lapse function depends only on t: = ^(t) 
[1] . Since the Hamiltonian constraint is derived from the 
variation with respect to the lapse function, in the pro- 
jectable version of the theory, the resultant constraint 
equation is not imposed locally at each point in space, 
but rather is an integration over the whole space. In the 
cosmological setting, the projectability condition results 
in an additional dust-like component in the Friedmann 
equation Q. 

The most generic form of the potential y-HL is given 

by 

Yhl = 2 A + .91 3?. 

+gslLA!ll + gc^V^^jkV'H"') , (2.6) 

where A is a cosmological constant, and % are the 
Ricci and scalar curvatures of the 3- metric gij, respec- 
tively, and giS (i = 1,...,9) are the dimensionless cou- 
pling constants. By a suitable rescaling of time we set 
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gi = —1. We also adopt the unit of = 1 (A^pl = 1) 
throughout the paper. 

Let us consider a Bianchi IX spacetime, which metric 
is written as 



a 

T 



where the invariant basis cj* is given by 

= — sinx^dx"'^ + sin.T^ cosx'^dx^ , 



= cosx'^d.T^ 



= cos a; dx 



sin X sin a; dx 



(2.7) 



(2.8) 



The basic equations describing the dynamics of Bianchi 
IX spacetime in HL gravity are now given by the foUow- 
ings: 



3(3A- 1) 



3(/3l+/3') + ^n«,/3±) + ^ 



3(3A- 1) 



da 



(2.10) 



0, (2.11) 



A typical scale of length of the universe is given by a, 
which reduces to the usual scale factor in the case of 
the FLRW universe. We shall call it a scale factor in 
Bianchi IX model as well. The traceless tensor (3ij mea- 
sures the anisotropy of the universe. The spacelike sec- 
tions of the Bianchi IX is isomorphic to a three-sphere 
S^, and a closed FLRW model is a special case of the 
above metric in the isotropic limit (Pij — > 0). 

For a vacuum spacetime, without loss of generality, we 
can assume that /3y is diagonalized as 

= diag (/3+ + V3I3^ , 13+ - %/3/?_ , -2/3+) . (2.9) 



/3.+3F/^, + g|1^0, (2.12) 



where H = d/a is the Hubble expansion parameter, i.e., 
the volume expansion rate is given hy K. — 3H. The 
constant C arises from theprojectability condition and 
it could be "dark matter" Q , but here we assume C = 
just for simplicity. 

The potential V, which depends on a as well as /3±, is 
defined by 



V{a,f3±) ^rHL = 14)(/3±) + V^i(/3±)a + y2(/3±)a' + 14(/3±)a" + ^a^ 



(2.13) 



where 



=.12/9+ 



4(ff5 + 356 + 57 - 859) cosh(12V3/3_) + ^(355 + 656 + 357 - 459) cosh(8V3/3_) 



4(1555 + 13.g6 + 1557 - 2459) cosh(4%/3/3_) + ^(655 + 356 + 657 - I259) 



^(355 + 556 + 357 - 459) cosh(10V3/3_) - ^(955 + 756 - 357 + 2O59) cosh(6%/3/3_ ) 
+ ^(355 +96- 3.97 + 12.99) cosh(2V3/3-) - ^(1555 + 1356 + I557 - 2459) cosh(8V3/3_) 
+ ^(355 + .96 - 357 + 1259) cosh(4V3/3_) + i(355 + 856 + 2757 - 7259) 

^(555 + 356 + 557 - 1259) cosh(6V3/3_) + i(355 + 56 - 357 + I2.99) cosh(2%/3/3-) 
+e-^2/3+ _ 1(15^5 + 1356 + 1557 - 2459) cosh(4y3/3-) - ^(855 + Tge - 357 + 2O59) 
-7(355 + 556 + 357 - 459)6-^^'"+ cosh(2V3/3-) - 4(55 + 356 + .97 - 859)6 



-6/3+ 



-2i/3+ 



54 

2 



(cosh(10V3/3_) - cosh(6V3/3_)) + e'^^+ cosh(8y3/3_) + cosh(4V3/3_)) 



ie-8/J+ _ e-i4^+ cosh(2V3/3_) + ig-^o/J, 



(2.14) 



(2.15) 
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{g2 + Sgs) cosh(8y3^_) - 4(^2 + 53) cosh(4\/3/3_) + 852 + 93 
-^[92 + 53)e'^+ [- cosh(6V3/3_) + cosli(2V3/3_)] + ^6-^"+ [(8.92 + 53) cosh(4V3/3_) +52+53 
7(ff2 +53)e-i°''+ cosh(2V3/3_) + ^-(52 + 393)e-'^^\ 



32 ~ 

- cosh(4V3/3_) + 1) + e-2'3+ cosh(2V3^_) - ^e^^'^^ 



(2.16) 
(2.17) 



Although the similar potential was found in |53|, we ex- 
tend it to the case without the detailed balance condition. 



III. LINEAR PERTURBATIONS OF THE FLRW 
UNIVERSE 

In this section, we shall analyze the present system 
with small anisotropics by linear perturbations of the 
FLRW universe. We discuss the stability of the oscil- 
lating FLRW universe and present a new type of non- 
singular solutions. 



A. Oscillating Closed FLRW Universe 

First we summarize the result of the closed FLRW 
spacctime, which metric is given by 



-df + 



We find the Friedmann equation as 
1 9 

-d" + W(a) = 0, 



r^dfl^ 



where 



«(«) 



1 



A , 
1 - -a^ 



9r 



9s 



3a2 3a4 



3A - 1 

The coefficients and gg are defined by 

5r := 6(53 + 852), 

5s := 12(955 + 356 + 57) ■ 



(3.1) 



(3.2) 



(3.8) 



(3.4) 



The conditions for an oscillating FLRW universe to 
exist were already given in which are summarized 
as follows: 



(a) A = 



3r > , - |l < < 



(3.5) 



(b) A = 8/^2 >o 



< 

< 9l'^''''-\9r) , 



(3.6) 



U{a) 

0,06 




(a) U{a) for A = 




(b) V{a) for A = 3/10 
U{a) 




(c) U{a) for A = -3/10 



FIG. 1: The potentials for the oscillating FLRW universe. 
We choose the coupling parameters as A = 1, = 6, and 
(a) 5s = -72/25, (b) = -72/25, and (c) = -66/25, 
respectively. 



(c) A = < 

5r > 0, 



5i-''''^"^(5r)<.9.<0, 



(8.7) 



where gr = 5r/^^, 5s = 9s/i'^, and 5. 

by 



~[i.i](±) j^-^-j jg defined 



9t'^'^^ i9r) 



9e2 



2-865^ + 2(1-65,) 



n ^3/2 



(8.8) 



with 6 ~ +1 being the sign of A. 
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We show the typical shapes of the potential 1/(a) in 
Fig. [1] for the coupling parameters which we use in 
our numerical analysis. For an oscillating universe, a 
is bounded in a finite range as Omin < a < Omax- 



B. Linear Stability of the FLRW Universe 

As we mentioned in Introduction, a non-singular 
FLRW universe such as an oscillating universe should be 
stable against anisotropic perturbations. Otherwise such 
a spacetime may not be realized in the history of the uni- 
verse. Hence, in this subsection we study stability of the 
FLRW universe against linear anisotropic perturbations. 

When the anisotropic part of metric vanishes, i.e., 
(3± = 0, Eq. (|2.10p reduces to the usual Friedmann equa- 
tion p.2p for a closed universe. For the stability analysis, 
we expand the potential V{a, /3±) around the FLRW uni- 
verse with /3± = to second order of f3± , so that 



(9.g5 + 3.g6 



16 
9 

~2 

954 {0+ 
3 



(955 

2 



56 



-.97) 
357 



459)(/3i + /3? 



(3.9) 
(3.10) 



^(53 + 352) + 2(53 - 352) {Pi + Pi) , (3.11) 



-Pi)] 



(3.12) 



The total potential V{a,l3±) is thus approximated by 

y(a,/3±) ~ Uo{a) + U2{a){l3l+pl) , (3.13) 

where 



C/o(a) 



U2{a) 



^(9.9a 



64 

9, 



356 
A 

— a 
64 



57 j 



^(53 + 352)0^ 



-2(^-95~56 - 357 - 4gg) 



954a 



-^(53 



3,92)0^ 



r 



(3.14) 



(3.15) 



In order to ensure that the FLRW universe is stable 
against small anisotropic perturbations, we impose the 
condition [/2(a) > 0. The sufficient condition for sta- 
bility is obtained if all coefficients in U2{a) are positive 
because a is positive, i.e., 



The necessary and sufficient conditions for stability are 
obtained by taking into account the dynamics of the 
background FLRW universe, i.e., the time evolution of 
a scale factor a. We will cover a wider range of the cou- 
pling parameters than the above. To show the explicit 
necessary and sufficient conditions, just for simplicity, 
we restrict our analysis to the case with f/4 = 0, i.e., the 
parity-conserved theory. In this case, U2 can be recast in 



U2{a) = - 



a^-l2g2 + ^] +48(56 +57 +.99) 



-5s - ( 1252 - I 



(3.18) 



The stability condition is such that 1/2(0) > for Va £ 
[oniin, amax] (the range in which the oscillation occurs). 
When A = 0, amin and amax arc given explicitly by 



V.9? + 12.9. 



5r + ^5? + 125. 



(3.19) 



We thus find that either of the following three conditions 
must be satisfied to ensure the stability: 



52 < 



^(35.-v^ 



5? + 12.9. 



56 + 57 + 59 > y^5s + ^ (2452 
X (sr - V.9? + 125s) , 



1.) 



(3.20) 



or 



52 > ^ (35,. + v/5? + 12.9e 



56 + 57 + 59 > ^5s 



288 



(2452 - 5r) 



X (fj, + v/5? + 125s) 



(3.21) 



or 



32-^5. 



< y^v/.9? + 125s 



56 



57 + 59 > i^5s + 7^ (3652 



48- 



432 



\2 
5r) 



(3.22) 



-955 ^ 

54 > 

53 



56 + 357 + 459 > , 



352 ^ - 652 > , 
6 

1 



52 S -^9r, 

36 

54 >0, 

1 

56 + 57 + 59 > -^9: 



(3.16) 



(3.17) 



The above inequalities give the necessary and sufficient 
conditions for a stable oscillating FLRW universe in the 
case of A = and 34 = 0. Since in more general cases 
with A 7^ and/or 34 7^ 0, the necessary and sufficient 
conditions could be obtained straightforwardly but would 
be much more involved, we dare not list the full condi- 
tions in the present paper. It may be sufficient to demon- 
strate that the stability range of coupling parameters ex- 
ists without any fine tuning. 
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C. Perturbation around a static universe 

Next we provide a simple and illustrative example in 
which even small anisotropies /3± can bring in a possibly 
interesting cosmological dynamics. 

Let us consider the case with + 3g2 > 0, (33 + 
3^2)^ + 4(9g5 + 3.96 + 57) = 0, and A = 0, so that we 
find a static FLRW universe with the constant scale fac- 
tor a = as '■= Vffs + 3.92 • We then add small anisotropic 
perturbations /3± to this static background. The basic 
equations governing the system arc given by 



= 

/3± + 



3A- 1 
32 dV 
3a6 dl3± 



fi4 

= 0, 



(3.23) 
(3.24) 



where 



y(a,/3±) ~ Uoia) + U2ia) {Pl + I3'i) . (3.25) 



Integrating Eqs. p.23p and p.24p . we find conserved 
anisotropic energies Ep^ defined by 



• 64 

Pl + —U2ias)l3l 



(3.26) 



Using the constants Ep^ we obtain the equation for the 
scale factor a as 



1 



(3.27) 



This equation gives an oscillating solution for a with the 
frequency 

(3.28) 



(3A-l)a|- 



Similarly, /3± also oscillate with the frequencies 
2 64C/2(as) 



3a| 



(3.29) 



The ratio of two frequencies tUa and uj/j^ is given by 



/8(3A- 1)^2(05) 
34 



Those frequencies give the typical values of the present 
oscillating system. They and their ratio are fixed only 
by the coupling parameters (gi's) because as and U2{as) 
are given by them. 

The above perturbative analysis around a static uni- 
verse shows that an oscillating solution newly appears in 
the presence of the small anisotropics f3±. As the argu- 
ment here is based on perturbations \l3±\ <ti 1, one may 
wonder whether or not there exist similar oscillating so- 
lutions with large anisotropies. We are going to perform 
numerical calculations to explore the anisotropic cosmo- 
logical dynamics arising from more general setups beyond 
perturbations. 



IV. ANISOTROPIC OSCILLATING UNIVERSE 
IN HORAVA-LIFSHITZ GRAVITY 

In the previous section, we have shown that an os- 
cillating FLRW solution in HL gravity is stable against 
small anisotropies, \f3±\ <C 1, for a wide range of the 
coupling parameters. We now proceed to investigate the 
rich variety of the dynamics of oscillating universes in the 
context of Bianchi IX spacetime, extending the analysis 
to the case with large anisotropies. 

If I3± are not small, the previous perturbative approach 
is no longer valid. To extend the analysis to include the 
case with large anisotropies, we employ a numerical ap- 
proach and solve the governing equations without any 
perturbative expansion. With this, we intend to uncover 
the rich variety of anisotropic cosmology and clarify the 
resultant fate of the universe. 

We have already given the basic equations for the 
Bianchi IX universe in HL gravity. It will be convenient 
to rewrite the equations as 



3(3A- 1) 



64 



8C 



a' + -Via,p±) + — 



H + iH = 



3(3A- 1) 
32 dV 



8 dV 
a^ da 



0, 



3C 



where 



(4.1) 
(4.2) 

(4.3) 
(4.4) 

(4.5) 



(o'ais) = diag((T+ -I- -\/3<T- , a+ — \/3a- , —2a+) is the shear 
tensor of a timelike normal vector perpendicular to the 
homogeneous three-space, and a is its magnitude. It may 
be convenient to introduce the dimensionlcss shear by 



S+ = — i and S = — 
H H 



(4.6) 



which measure the relative anisotropies to the expansion 
rate H. We also introduce the phase variable (p defined 

by 



arctan 



0-+ 



(4.7) 



which parameterizes the direction of the anisotropic ex- 
pansion. 

We have the five first-order evolution equations for H, 
l3± and a±, i.e., Eqs. (P3|)and supplemented 

with one constraint (|4.ip . We have to set up the initial 
values for five of the six variables, a,i?, /3±,E^ and cp. 
The other one is fixed by the constraint equation. Since 
we are interested in how the cosmological dynamics is 
altered by the introduction of anisotropies, we start with 
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the isotropic oscillating universe by setting /3±^o to vanish 
and Oq to be a local minimum of Zi(a), with arbitrary 
shear (Eq and (po) at the initial moment. So we shall 
give the initial data for ao, /3±.o, Eq'^ and (pa, and then 
determine Ho (or do) by the constraint (|4.ip . 

Without any loss of generality, we can analyze only the 
range of < (/Jq < 'I'/S because of the discrete symmetry 
modulo 7r/3 of the potential V{a,f3±). Since U{aM) is 
negative (therefore T^(aAf,/3± = 0) is positive) for the 
oscillating FLRW universe, the possible range of initial 
shear Sg is limited from Eq. (|4.2p as 

< < Sg,„a.) ■■= ^^^^ ■ (4.8) 



As we discussed in the previous section, the FLRW 
universe is stable against small anisotropic perturbations 
if U2{a) defined by Eq. p.l8|) is positive. However one 
may suspect that it becomes unstable when anisotropy 
is large. For stability against large anisotropy, we have 
one natural indicator, which is the potential V{a^ P±)- If 
the potential is unbounded from below for large |/3±|, we 
expect that if initially large | /3± | will diverge in time and 
the universe evolves into a singularity. 

From Eq. (|2.14p , we find that the potential is bounded 
from below, if and only if one of the following conditions 
is satisfied: 



(i) 59 > g (97 + 3.96 +35), (4.9) 

(ii) 54 > 0, bgj + 7.96 + 5.95 > , and .99 = ^{gr + 3^6 + 55) with gj + (557 + 7^6 + ^gs f + , (4.10) 

(iii) 3.93 + 32 > 0, 36 > ,.94 = 0, .97 = -\{'^9^ ~ ^9^) , and 39 = ^36 with (333 + 32)^ + 36 7^ , (4+1) 

(iv) .93 > 0, 35 < , .92 = -333, .97 -35 , and 54 56 = 39 = . (4.12) 

I 



We then classify the potential V into four types: SS, 
US, SU, and UU, where the first S (stable) or U (un- 
stable) denotes the stability against small perturbations 
around FLRW spacetime, while the second S (stable) or 
U (unstable) corresponds to the stability against large 
anisotropies. Since a cosmological constant will also 
change the fate of the universe, we shall discuss eight 
types of cosmological models; Models I-SS, I-SU, I-US, 
LUU, n-SS, n-SU, n-US, and n-UU, depending on the 
sign of A (I for A < and II for A > 0) and the potential 
types. 

Note that a singularity, of course, may appear even 
for small anisotropies because the Bianchi IX spacetime 
includes a closed FLRW model. 

We have performed numerical calculations for all pos- 
sible models and various initial data. Now we show our 
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TABLE I: The values of nontrivial coupling parameters gi's 
and a cosmological constant A, which are used in our numer- 
ical analysis. We also choose A = 1, 512 = fl4 = fl7 = 0. The 
types (SS and SU) of the potential V are described in the 
text. 



numerical results for each model. We find several types 
of fates of the universe depending on the magnitude of 
anisotropy, which we shall describe one by one. 

In Table U we list up the values of parameters for which 
we present the figures in this paper. 

A. Model I-SS (A < and Type-SS potential) 

First we discuss Model I-SS, in which A < and the 
potential V is the SS type. Since the universe is closed, 
there are two fates: an eternal oscillation or a big crunch. 
Depending on the strength of initial anisotropy, we find 
the following three types of histories of the universe. 

(A) Anisotropic oscillation : (small anisotropy) 

Since we have an oscillating universe for the FLRW 
spacetime, we find an eternally oscillating non- 
singular solution if Eq is sufficiently small (Eq < 
^o(max)/>^)- ^ typical example is given in Fig. [21 
The scale factor a is regularly oscillating with time 
just as the FLRW solution with the same initial 
scale factor ao, which is shown by the dotted red 
curves as a reference. 

This oscillating solution shows only small deviation 
from the isotropic FLRW universe. The scale factor 
a (and then the volume) oscillates very regularly. 
Its orbit in the phase space shows an ellipse (a cross 
section of a torus) (see Fig. Hlb)). The radius 
is slightly larger than that of the FLRW universe 
because of the existence of shear (see Eq. (14.11) ) 
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(a) the time evolution of scale factor a 
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(b) the phase space of (a, d) 

FIG. 2: The time evolution of non-singular oscillating uni- 
verse. The solid blue line and dashed red line represent 
Bianchi IX universe and FLRW universe with the same initial 
ao, respectively. We show the scale factor a in (a) and the 
orbit of a in the phase space in (b), respectively. We set the 
initial values as ao = 2^/6/5, ho = 0.2282, Eg = 0.6000 and 
ifo = 2tv/9. 



The orbit of the anisotropy is depicted in 

Fig. ^a). The anisotropic variables f3± arc trapped 
around the origin of /3_)-space by the poten- 
tial wall. It looks complicated but definitely pe- 
riodic. The shear is also regularly oscillating as 
shown in Fig. [St^b), but the oscillation period is 
much shorter than that of the scale factor a. The 
oscillation amplitude of the shear a is then modu- 




(a) the orbit of 




(b) the evolution of the shear 



FIG. 3: The orbit of anisotropy (/3+, /3_) and the time evolu- 
tion of the shear of the oscillating universe given in Fig. 



lated by a-oscillation. We can estimate those oscil- 
lation frequencies from the result in mil CI Using 
the coupling parameters of the present model (see 
Table [I]), we find cJq = 2 and = 2-\/T0 from 
Eqs. p.28p and p.29p . which are almost the same 
as the frequencies in Figs[2Ja) and[3l^b). 

The resultant universe is regarded as an isotropic 
spacetime with small anisotropic perturbations. 
The anisotropic oscillation may continue eternally. 



(B) Big crunch after oscillations : (large anisotropy) 

If En is large as E^, >/3 < En < En, an ini- 
tially oscillating universe eventually collapses into a 
big crunch (a = 0) after many oscillations because 
of increase of the anisotropy. A typical example of 
such a singular universe is shown in Fig. 21 
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(b) phase space (a, d) 



FIG. 4: The time evolution of the unstable oscillating Bianchi 
IX universe. The dashed red line shows the FLRW spacetime. 
We set the initial values as ao — 2y/Q/5^ do = 0.5455, Eq = 
2.5800 and ipa = 27r/9. 

The oscillation period is almost the same as that in 
Fig. lUJa). As shown in Fig. S] (b), the orbit of the 
scale factor a in the phase space is initially almost 
an ellipse (a cross section of a torus), but its "ra- 
dius" gradually increases because of the increasing 
anisotropy, and a finally evolves into a big crunch 
singularity ((a, d) = (0, — oo)). 

The behaviour of anisotropy is shown in Fig. [5ja) , 
which shows that the orbit of (/3+ , /3- ) is trapped 
and reflected many times by the potential wall. The 
shear a is initially oscillating and eventually di- 
verges at a big crunch as shown in Fig.[SJb). Before 
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(a) the orbit of (/?+,/?-) 




(b) the evolution of shear cr^ 



FIG. 5: The time evolution of the shear of the unstable oscil- 
lating universe given in Fig. 2] 

this divergence, we can see the increase of cr^ , which 
leads the leave from the oscillating phase. Then the 
universe eventually evolves into a singularity with 
finite values of (3±. Note that the relative shear S 
is finite at the end, which means that the shear is 
not responsible for the singularity. 

In Fig. ini we show the curvature invariant 

^y^'J =3i^2 + 2a^ (4.13) 

which really diverges at a big crunch singularity. 
The universe evolves into a big crunch after many 
oscillations. 
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FIG. 6: The time evolution of the extrinsic curvature square 
5C %ii of the unstable oscillating universe given in Fig. |4l 
For the Bianchi IX universe (solid blue line), 5C oscil- 
lates in the beginning, but it eventually diverges at the end 
of the evolution, while it just oscillates periodically for the 
FLRW universe. 

(C) From big bang to big crunch : 

(near maximally large anisotropy) 

Another type of singular solution is found for the 
extremely large initial anisotropy (Eq ~ ^o(max))- 
In the case of the closed FLRW universe in GR, the 



spacetime starts from a big bang and ends up with a 
big crunch. Bianchi IX universe in GR also evolves 
from zero volume (a big bang) to zero volume (a big 
crunch) through a finite maximum volume. Hence 
even for the case with the oscillating FLRW uni- 
verse, if we add a sufficiently large anisotropy, we 
may expect such a non-oscillating simple evolution. 

We show one example. As shown in Fig. [71 the 
scale factor evolves from an initial finite value oq to 
a singularity (a = 0), which is called a big crunch. 
If we calculate the time reversal one from the same 
initial data, we will find a singularity (a = 0), which 
is called a big bang. As a result, this universe starts 
from a big bang and ends up with a big crunch. 
There is no oscillation in the evolution of a just as 
the closed FLRW universe in GR. 



a 




FIG. 7: The time evolution of the scale factor a of the Bianchi 
IX universe with the potential given in Fig. 12.61 is shown 
by the solid blue line. The dashed red line represents the 
oscillating FLRW universe. We set set the initial values as 
ao = 2^/6/5, ao = 2.0412, Eg = 0.990E(^(„,^,,) (= 2.9700) and 
ipo = n/6. 



As for the anisotropy, as shown in Fig. |8l the orbit 
of f3± is oscillating around the origin and refiecting 
at the potential wall until formation of a singularity. 
f]± is finite even at a big crunch. 

The shear a is also oscillating, but the frequency 
is not so regular compared with the previous two 
cases (A) and (B). 

The shear diverges at a big crunch (a = 0), which 
is really singular because the extrinsic curvature 
square ^ij also diverges there as shown in Fig. 
[9l The behaviour of ^ '-^ ^ ij is very similar to that 
of the shear square. However the relative shear 
E does not diverge at a big crunch, which means 
that the singularity is similar to that of the FLRW 
universe. The shear does not dominate in the dy- 
namics (Compare it with next example (C)'). 

In the case of a negative cosmological constant 
(A < 0), we also find the similar behaviour of the 
universe, although there is a quantitative differ- 
ence. 
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(b) shear square 



FIG. 8: (a) The orbit of anisotropy (/3+, /3_) and (b) the time 
evolution of the shear square of the universe given in Fig. 
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-5 




FIG. 9: The time evolution of the extrinsic curvature square 
?C %ii for the solution shown in Fig. [T] The solid blue line 
and dashed red line represent Bianchi IX universe and FLRW 
universe, respectively. It diverges at the end of evolution, 
which is a singularity. 

B. Model I-SU (A < and Type-SU potential) 

When the potential V is unbounded from be- 
low, the universe may be unstable against large 
anisotropic perturbations (see Fig. [TU)) . We also 
find three types of histories of the universe as (A), 
(B), and (C) in Model I-SS. The difference from 
Model I-SS appears when anisotropy gets large. 
That is, the anisotropy /3± diverges when a sin- 
gularity appears. We show one example with near 
maximally large initial anisotropy. 

(C)' From big bang to big crunch: 

(near maximally large anisotropy) 

In this case, just as the history (C) of Model I-SS, 
the Bianchi IX universe expands from zero volume 
to zero volume without oscillation (Fig. [TT|) . 




FIG. 10: The unstable potential V{a,j3±) against large 
anisotropic perturbations is shown for a = 1. 



a 




1 2 3 4 5 6 



FIG. 11: The time evolution of the scale factor a of the 
Bianchi IX universe is shown by the solid blue line. The 
dashed red line represents the oscillating FLRW universe. 
We set set the initial values as ao = 2^/6/5, ho = 6.4549, 
Eg = 2.9970 and ipo = vr/lS. 



The difference appears in the behaviour of 
which diverges at a big crunch. The orbit of I3± 
initially oscillates around the origin and reflects on 
the potential wall, but it eventually evolves to in- 
finity over the potential hill because the potential 
is not bounded from below (see Fig. [T^a)). We 
also show the time evolution of the shear cr^ , which 
diverges at the end of the universe. The extrin- 
sic curvature square Ki 'K.ij also diverges there, 
which means it is really a singularity. 

This singularity is different from one appeared in 
the history (B) or (C). To show it, we depict the 
time evolution of the relative shear S, which di- 
verges at a big crunch. It means that the shear 
becomes dominant at the end. The increase of the 
anisotropic shear is responsible for the formation 
of a singularity. Note that S diverges also in the 
middle of the evolution but its divergence appears 
because of = 0. 

At a big bang, which appears in the time reversal 
one, we suspect that the shear diverges but /3± is 
finite just as the beginning of Bianchi IX universe 
in GR. 
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(a) anisotropy (/3+,/3_) 
log(a2) 




(b) shear square 



FIG. 12: (a) The orbit of anisotropy (/3+,/3-) and (b) the 
time evolution of the shear cr^ of the universe given in Fig. 
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FIG. 13: The time evolution of the relative shear square E 
of the universe shown in Fig. m 



C. Model II-SS (A > and Type-SS potential) 

Next we discuss the case of A > 0. In this case, we find 
another fate of the universe, which is an exponentially 
expanding universe by a positive cosmological constant. 

If the effect of anisotropy is smaller than the contri- 
bution of a cosmological constant, the universe will be 
isotropized. The asymptotic equation for a is given by 



a 6\ 



2A 



1 



a. 



(4.14) 



if we neglect the anisotropic terms in Eq. (|2.1ip . find- 
ing an exponentially expanding FLRW spacctime, i.e., 
de Sitter spacetime; 



(4.15) 



Along with the cosmic expansion the potential is flat- 
tened to be negligible in the evolution equation for the 
anisotropy because a increases rapidly, so that we find 



(4.17) 



from Eq. p.24|) . This implies that ~^ Q asymptot- 
ically due to the Hubble friction. Since the shear de- 
creases to zero, the universe becomes locally isotropic, 
i.e., locally de Sitter spacetime. Note that this does not 
mean that it is globally de Sitter spacetime because the 
asymptotic values of /3± do not vanish. However the 
spacctime is exponentially expanding, and the observ- 
able region such as a horizon scale becomes effectively 
isotropic. Hence we shall still call this asymptotic space- 
time de Sitter universe. 

As a result, we have three fates of the universe in the 
case with A > 0: oscillating universe, a big crunch, and 
dc Sitter expanding universe. We then find five types of 
histories of the universe: two new types with asymptoti- 
cally dc Sitter universe in addition to the previous three 
types (A), (B) and (C) discussed in IIV Al 

Two new types are the similar to the histories (B) and 
(C), but different from those in their final states, i.e., (D) 
de Sitter expansion after oscillation, and (E) de Sitter 
expansion from a big bang. 

For small anisotropy (E§ < Sp^j^^^^/S), wc find the os- 
cillating universe (A) . While for the large anisotropy 
(So(max)/3 <^l< So(„,ax)). t^e universe is either Type 
(B) or Type (D). When the initial anisotropy is extremely 
large enough, i.e., <So(„,^^), Type (C) or Type (E) 
is obtained. We shall describe new types (D) and (E) 
below. 



(D) de Sitter expansion after oscillation : 
(large anisotropy) 

We show one example in Fig. [HI 




FIG. 14: The time evolution of the scale factor a of the uni- 
verse with A > 0. The solid blue line and dashed red line rep- 
resent Bianchi IX universe and FLRW universe, respectively. 
We choose the coupling constants as gs = 1, = —3/100, 
gg — 1/100,32 = (/4 = (?6 = 57 = 0, A = 3/10, and A = 1, and 



with 



set the initial values as ag 
and (fio — 7r/9. 



1.0758, ho = 0.6125, = 2.7000 



2A 



3(3A- 1) 



(4.16) 



If S2 is as large as S^^^^^^/S < < E^^^^^^, 
an initially oscillating universe eventually evolves 
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(a) 




(b) /34 




(c) /3_ 



FIG. 15: The time evolution of the shear square and 
anisotropies /3± of the universe shown in Fig. 1141 The shear 
oscillates initially, but the universe suddenly leaves the oscil- 
lating phase to de Sitter expanding phase, which drops the 
shear to zero rapidly. The initially oscillating P± finally set- 
tles to finite values after small bump. 



into an exponentially expanding dc Sitter universe 
because of a cosmological constant. 

The initially oscillating universe leaves the oscilla- 
tion phase when the anisotropy increases beyond 
some critical value. We also show the evolution of 
the shear a in FigfTSl We find that it is oscillat- 
ing regularly for two-third of the whole period, but 
eventually increases. Then the universe leaves the 
oscillating phase and evolves into de Sitter phase. 
Because of rapid expansion of the universe, the 
shear vanishes soon [6l|. 

We also show the time evolution of the anisotropy 
(/3±) in Fig. [T5l b).fc). The initially oscillating 
anisotropy increases as a burst and then decreases 
to a small finite constant. 

How the universe choose its fate ((B) or (D)) is 
as follows: If the spacetime is expanding when it 
leaves from the oscillating phase, it evolves into de 
Sitter phase (D), while if it is contracting, it col- 



lapses to a big crunch (B). 



(E) de Sitter expansion from big bang : 
(near maximally large anisotropy) 

If Eq is close to the maximum value (2^^.^^^^^^^^ = 3) 
and the universe is initially expanding (ag > 0), the 
spacetime evolves into de Sitter phase without os- 
cillation. The large anisotropy makes a jump from 
the oscillating phase to de Sitter phase in the begin- 
ning. The anisotropy drops quite rapidly because 
of the exponential expansion as shown in Fig. 1171 




FIG. 16: The time evolution of the scale factor a of the uni- 
verse with A > and Type-SU potential. The universe starts 
from a big bang and evolves into de Sitter spacetime. The 
dashed red line represents the oscillating FLRW universe as 
reference. 

On the other hand, if the universe is initially con- 
tracting (do < 0), the spacetime is classified into 
Type (C), i.e., from a big bang, which appears in 
the time reversal one, to a big crunch without os- 
cillation. 



log(cr 




FIG. 17: The time evolution of the shear of the universe 
shown in Fig. 1161 Initially oscillating shears drops to zero 
after de Sitter expansion starts. 



D. Model II-SU (A > and Type-SU potential) 

In this case, we also find the similar histories of 
the universe to Types (A), (B), (C), (D) and (E), 
depending on initial anisotropies. The differences 
between Models 11-SS and 11-SU are qualitatively 
the same as those between Models I-SS and I-SU. 
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Only one difference from Model I is that there exists 
de Sitter phase as the fate of the universe because 
of a positive cosmological constant. 



E. Models with the unstable potential against 
small perturbations around FLRW spacetime 

In the previous four subsections, we discuss the 
cosmological models with the stable potential 
against small perturbations around FLRW space- 
time. When the potential is unstable against small 
perturbations around FLRW spacetime, we also 

I 

F. Dependence of anisotropy on the date of the universe 

In Table HIl we summarize the fate of the universe. We assume the coupling parameters by which there exists 
an oscillating FLRW universe. For Models I-SS and II-SS, we find an oscillating FLRW universe with anisotropy in 
the case of small initial anisotropy. When we increase the strength of anisotropy, the spacetime leaves the initially 
oscillating phase and eventually evolves into a singularity or de Sitter spacetime. If the initial anisotropy is near the 
maximum value, the oscillating phase disappears and a simply expanding and contracting universe is found just as a 
closed universe in GR for A < 0. When A > 0, an initially expanding universe evolves into de Sitter spacetime, while 
an initially contracting universe evolves into a big crunch. 



Model 


cosmological 


potential 






constant 


V 


small 


t> 


large 


[> 


near maximally large 


I-SS 




SS 


(A) OSC 


> 


(B) OSC-^SINGl 


O 


(C) SINGl 


I-US 


A < 


US 


(A) OSC 


> 


(B) OSC^SINGl 


O 


(C) SINGl 


I-SU 




su 


(A) OSC 




> 




(C)' SING2 


I-UU 




uu 






(C)' SING2 






II-SS 




SS 


(A) OSC 


> 


(B)-(D) OSC'^deS/SINGl 


> 


(E) deS [(C) SINGl] 


II-US 


A > 


us 


(A) OSC 


> 


(B)-(D) OSC'-^deS/SINGl 




(E) deS [(C) SINGl] 


II-SU 




su 


(A) OSC 




> 




(E)' deS [(C)' SING2] 


II-UU 




uu 






(E)' deS [(C)' SING2] 







TABLE II: Classification of Bianchi IX cosmological models by a cosmological constant A and types of the potential. OSC, 
SINGl, SING2 and de S represent the anisotropic oscillation, two types of big crunch singularities (one with finite anisotropy 
and the other with infinite anisotropy), and de Sitter spacetime, respectively. deS/SINGl means that the spacetime evolves 
either de Sitter phase or a big crunch with finite anisotropy. deS [SINGl or SING2] denotes that the fate is either de Sitter 
universe if the universe is initially expanding or a big crunch if contracting. 



find qualitatively similar results. The main differ- 
ence is that oscillations around the FLRW space- 
time never happen. Even if the universe starts from 
near FLRW spacetime, it evolves into spacetime 
with large anisotropy because the FLRW space- 
time is unstable. As a result, in the case of Type- 
UU potential, the universe collapses to a singularity 
for Model I-UU. No oscillating phase is found. If 
A > 0, i.e., for Model II-UU, some universe col- 
lapses to a singularity without oscillations, and the 
other one evolves into de Sitter expanding universe, 
depending on initial conditions. 



For Models TUS and II-US, the histories of the uni- 
verses are similar to those in Models I-SS and II-US, 
respectively, although deviation from isotropy becomes 
large even for initially small anisotropy. 

In the cases of Models I-SU and II-SU, unless an initial 
anisotropy is small, an initially expanding universe turns 
to contract and collapses into a big crunch for A < 0, 
while it evolves into de Sitter spacetime for A > 0. This 
singularity at a big crunch is different from that in Mod- 
els TSS,-US, and ITSS, US. The anisotropy /3± for the 
present model diverges, while that for the other cases is 



finite even at a singularity. 

Models TUU and II-UU are not so interesting. There 
is no oscillating phase. If A < 0, the spacetime is simply 
a spacetime evolving from a big bang to a big crunch. 
For A > 0, the initially expanding universe evolves into 
de Sitter spacetime, while the contracting one collapses 
into a singularity. 

One may wonder whether the initial anisotropy classi- 
fies the fate of the universe. Are there any critical values 
of the initial anisotropics for their transitions in Table |TT] 
? We understand naively that such a transition occurs 
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as the anisotropy increases. 

We find there exists a critical value for the transition 
(B) to (C) or (D) to (E), which is about (0.94- 0.98) x 
^o(max)- transition from (A) to (B) (or (D)), it is 

not so clear whether there exists a critical value of Sq or 
not. If Eq is sufhcicntly small, we find the history (A), 
while when it is large, we find the history (B) or (D). 
However, because we analyze the system numerically, we 
are not sure whether the model with small anisotropy 
oscillates forever or will turn to collapse long after. If 
the latter case is true, the case with small anisotropy is 
classified into the history (B) or (D). So there may not 
be exactly the history (A) except for the exact FLRW 
spacetimc. 

More interesting fact is found in the history (B)-(D) for 
A > 0. In order to study how the fate of the universe de- 
pends on the initial data, we solve the basic equations in 
Model II-SS assuming various initial values of anisotropy 
(Eq and i^o)- We summarize the results in Fig. [TH] As 
we see from Fig. [T8j the spacetime with do > start- 
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FIG. 18: The fate of Type II-SS Bianchi IX universe in terms 
of initial anisotropy Eq and (po- We judge the fate of the uni- 
verse at t = lOOtpL {tpL'- the Planck time). The histories (A), 
(B), (C), (D) and (E) are represented by a filled green circle, 
yellow triangle, red cross, purple diamond and blue square, 
respectively. The empty yellow triangle and purple diamond 
are classified to the histories (B) and (D), respectively, but 
those oscillating periods are longer than lOOtpL . We have set 
ao = 1.0758. do is fixed by the constraint equation [(a) do > 
and (b) do < 0]. 



ing from near maximum anisotropy evolves into de Sitter 
universe (the history (E)), while it collapses into a sin- 
gularity if do < (the history (C)). For small Eq, we 
find an oscillating universe (the history (A)). If Eq is 
between the above two cases, however, the fate of the 
universe is not so simple. The history of such a universe 
is classified either (B) OSC-^SINGl or (D) OSC^deS. 
However such a history does not shift monotonically from 
(B) to (D) as the initial anisotropy increases. How the 
universe choose its fate ((B) or (D)) as follows: If the 
spacetime is expanding when it leaves from the oscillat- 
ing phase, it evolves into de Sitter phase (D), while if it 
is contracting, it collapses to a big crunch (B). As a re- 
sult, the fate of the universe is sensitively dependent on 
initial conditions. If one changes the initial conditions, 
the fate changes drastically. It is because the present 
system is iioii-integrable. Such a property is found in a 
dynamical system with chaos. Since our model is Bianchi 
IX, which shows chaotic behaviour near singularity in 
GR [55| - |57| . we understand why we find such a compli- 
cated basin structure of the fate in Fig. [THl which can be 
fractal [57 1 . 



V. SUMMARY AND REMARKS 

We have explored a singularity avoidance in a vacuum 
Bianchi IX universe in HL gravity. We have studied an 
oscillating cosmological solution with anisotropy. In the 
case of small anisotropy (|/3±| ^ 1), we find an analytical 
solution and show the stability condition of the FLRW 
universe against anisotropic perturbations. We have also 
solved the basic equations numerically and discussed the 
possible history of the universe. We classify our models 
into eight types I-SS,-SU,-US,-UU, and II-SS, -SU,-US,- 
UU, depending on the sign of a cosmological constant 
and the types of the potential V. Wc find five types of 
the histories of the universe: (A) an oscillating universe 
with anisotropy, (B) a big crunch after oscillations, (C) 
from a big bang to a big crunch, (D) de Sitter expansion 
after oscillations and (E) from a big bang to de Sitter 
expansion, as summarized in Table HIl 

The stable oscillating universe (A) is found if initial 
anisotropy is small in the case that the coupling param- 
eters {gi 's) satisfy the stability condition. When initial 
anisotropy is large, the oscillating universe evolves into a 
singular big crunch (B) for A < 0. In the case of A > 0, 
if the initial anisotropy is large but not close to the max- 
imum value, we find two histories (B) and (D). Which 
history is realized does not depend monotonically on the 
initial shear Eq, but the present system shows sensitive 
dependence of initial conditions just as one of the typical 
properties of chaos. The anisotropic bounce universe is 
also obtained for the model satisfied the stability condi- 
tion if initial anisotropy is small. 

Since we adopt the unit of AfpL = 1, the oscillation pe- 
riod and oscillation amplitude are the Planck scale, unless 
the coupling constants arc unnaturally large. Hence in 
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order to obtain a macroscopic universe, we need a posi- 
tive cosmological constant A > 0, which provides us a de 
Sitter expanding phase. 

In a more reahstic situation, this cosmological constant 
should be replaced by a potential of an inflaton scalar 
field (j). Reheating after inflation may give an initial state 
of a macroscopic big bang universe. 

When we include a scalar field, however, we have to 
take into account modification of a scalar field action in 
the UV limit similar to HL gravity action. The action 
may be given by 



where 



dtd^xNy/g 



2N^ 



^0<t>-V^{dp) , (5.1) 



+ CiA 



(5.2) 



M4 AP 

with M being a typical mass scale and Ci's {i — 1 ~ 3) 
being dimensionlcss constants. For this action with HL 
gravity action (|2.ip . the basic equations for a Bianchi IX 
cosmological model are given by 



HA 



3(3A- 1) 
8C 1 -2 
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3H 



V4,{<l>) , 
8 dV 



3(3A- 1) [a^ da 



3C 



(5.3) 
,(5.4) 
(5.5) 
(5.6) 



Although the action (|5.ip contains higher spatial deriva- 
tives, there exists no difference from the conventional 
canonical kinetic term of a scalar field for a homogeneous 
spacetime. We expect a usual inflationary scenario once 
de Sitter exponential expansion starts. There exists re- 
heating after slow-roll inflation, finding a big bang uni- 
verse. 

If we have an oscillating phase before inflation, we may 
expect one interesting effect, which is modification of pri- 
mordial perturbations. We should stress that a classical 
transition from an oscillating phase to an inflationary 
stage never happens in the FLRW model (See our discus- 
sion in [47t for quantum transition). So anisotropy may 
be important in the pre-inflationary oscillating phase. As 
a result, we may find large no n- Gaussian density per- 
turbations. To confirm our scenario, we should explore 
the dynamics of a scalar field in the prc-inflationary era 
because in an inflationary model may be fluctuated 
by the oscillating scale factor in the pre-inflationary era 
while A is a constant in our present analysis. 

Another interesting possibility is anisotropic inflation, 
which is discussed in the model with higher curvature 
terms [6^ or with a vector field [63j. It may leave distin- 
guishable imprints on the primordial perturbations. 
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Appendix A: Anisotropic bounce universe in 
Hofava-Lifshitz gravity 

In the text, we have considered oscillating universes 
as a possible way to avoid an initial singularity. In this 
appendix, we study another way to singularity avoidance, 
i.e., a bounce solution [l^l- For a closed FLRW universe, 
a bounce solution exists only only in the case of A > 0. 
As a result, a bounce solution with anisotropy is also 
found only for the case of A > 0. 

We classify bounce solutions into two classes according 
to the shape of the potential U{a). The first class is 
such that U{0) = — oo, which we call Type A, and the 
other is such that U{0) = +oo, which we call Type B. 
We show two typical potentials in Fig. [121 For Type- 




(a) Type A 



U{a) 








/U.5 2.0 


2.5 \ 3.1 




^max 





(b) Type B 



FIG. 19: Typical shapes of potential W(a). (a) Type A: We 
choose the coupling constants as g,- = —6, gs ~ 27/25 and 
A = 3/10. The bounce point in the isotopic FLRW model is 
given by ar = 3.4183. (b) Type B: We choose the coupling 
constants as g^ = 6, (?s = —81/25 and A — 3/10. The bounce 
point, the maximum and minimum oscillating radii in the 
isotopic FLRW model are given by ar = 2.7452, flmax = 
1.2337, and amin — 0.9702, respectively. 

A potential, in addition to a bounce solution, wc find a 
FLRW spacetime starting from a big bang and collapsing 
to a big crunch, while for Typc-B potential, we have an 
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oscillating FLRW solution as well as a bounce spacetime. 
This difference gives rise to different fates of anisotropic 
Bianchi IX universe as we will show later. 

When initial anisotropy is small, it is easy to find 
anisotropic bounce solutions because anisotropy can be 
treated as perturbations around the isotropic closed 
FLRW universe. In order for this type of stable solu- 
tions to exist, the condition U2{a) > (Va > ax) must 
be satisfied. The sufficient condition for stability is the 
same as Eq. p.l7|) . 

If the stability condition U2{a) > is not satisfied, 
we do not usually obtain bounce solutions. Only for ex- 
tremely small initial shear (Sg < 10~^), a bounce solu- 
tion can be found even if the condition is not satisfied. 
It is because the the bounce occurs before the unstable 
mode grows enough. 

In Table IIIIl we list up the values of parameters for 
which we present the figures and tables here. 
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Figs. 18(a), 19, Table IV 
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Figs. 18(b), 20, Table V 



TABLE III: The values of nontrivial coupling parameters Qi's 
and a cosmological constant A, which are used in our numer- 
ical analysis. We also choose X = 1, g2 = = gr = 0. The 
types (A and B) of the potential U are described in the text. 



1. Type- A bounce universe 

First we show the results for Type-A potential in 
Fig. [201 where we find two typical evolutions of the uni- 
verse. 

log(a) 



















1— ' 


2 


3 


4 5 



FIG. 20: The time evolutions of the scale factor a of the 
universe with Type-A potential. The solid blue, and dashed 
red lines denote de Sitter universe via a bounce, and a col- 
lapsing universe, respectively. We have set ao — 3. 7601, and 
ipo ~ 57r/18, and the shear square Eq — 0.6900 (solid blue) 
and Eq — 1.9500 (dash-dotted red), do is determined by 
the constraint equations as do = —0.5968 (solid blue) and 
do = -0.8852 (dash-dotted red). 



We assume the coupling parameters which guarantee 
the existence of the FLRW bounce universe. We then 
include anisotropy and study how anisotropy changes the 
fate of the universe. We find the following results: 

• If the shear is small, we still have a regular bounce 
solution (the solid blue curve). 

• When the shear becomes large, the initially con- 
tacting universe just collapses to a singularity (the 
dash-dotted red line). No singularity avoidance is 
obtained. 

Hence, there seems to exist a critical value of 1]^^^^,^^^, 
beyond which the spacetime collapses to a singularity. 
The critical values Sp^^^.;^^ for Type-A potential are listed 
in Tabic IIVI which is strongly dependent on the initial 
scale factor oq. However, the corresponding values of 
^o(crit) ^^'^ ^^^^ ™uch different. Hence we may conclude 
that the critical value is determined by the absolute value 
of the shear but not by the relative value to the Hubble 
parameter. 



ao 


shear 


<^0 


ri 7r TT 7r 27r 57r 7r 
" 18 9 6 9 18 3 


l.laT 


'^O(crit) 
^2 

"O(crit) 


1.539 0.963 0.861 0.786 0.738 0.708 0.699 
0.426 0.285 0.264 0.246 0.237 0.231 0.228 


1.5aT 


y.2 

^O(crit) 
^2 

''O(crit) 


0.225 0.114 0.102 0.093 0.087 0.081 0.081 
0.207 0.144 0.135 0.129 0.126 0.120 0.120 









TABLE IV: The critical relative shear square Hqi^^j.^^-^ and 
the corresponding absolute shear square o-qi^^^^^-j for Type-B 
bounce universe, for Type-A bounce universe. If the initial 
shear exceeds this critical value, the universe evolves into sin- 
gularity. 



2. Type-B bounce universe 

In this case, we find the following three evolutionary 
histories of the universe: 

• A simple bounce solution just as Type A 

This is possible if the deviation from the "back- 
ground" isotropic universe is sufficiently small. 

• A big crunch solution just as Type A 

The universe collapses into a singularity if Sg is 
larger than some critical value Sgj.^j.j^^ . 

• A bounce solution after some oscillations 

Wc also find an oscillating phase before a bounce, 
for which Sq is very close to ^Qr^^■^^) initially. This 
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log(a) 




FIG. 21: The time evolutions of the scale factor a of the 
universe with Type-B potential. The solid blue, dash-dotted 
red, and dotted green lines denote de Sitter universe via a 
bounce, a collapsing universe, and de Sitter universe after 
several oscillations, respectively. We have set ao = 3.0197, 
ifio = it/9, and the shear square Eg = 0.9000 (solid blue). 
Eg = 1.1700 (dash-dotted red), and Eg = 1.1409 (dotted 
green), ao is determined by the constraint equations as ao = 
-0.4109 (solid blue), do = -0.4402 (dash-dotted red) and 
do = —0.4367 (dotted green). 

type of solution requires a fine-tuning to some de- 
gree, and hence in this sense the solutions are not 
generic. 



The typical evolution of Type-B universe is presented 
in Fig. [2TJ The critical values for Type-B bounce 

universes are listed in Table fVl As Type-A bounce solu- 
tion, the critical value of crQ(.^j,jj-) does not strongly depend 
on the initial scale factor ao, but S^^^^j^^ does. 

Near the critical anisotropy, we suspect that which fate 
of the universe is realized may depend sensitively on ini- 
tial conditions. 



ao 


shear 




ri 7r 7r 7r 27r Stt tt 
" 18 9 6 9 18 3 


l.laT 


^O(crit) 
^2 

'-^a(crit) 


1.254 1.188 1.143 1.086 1.065 1.062 1.137 
0.288 0.276 0.267 0.255 0.252 0.252 0.264 


1.5aT 


y2 

^O(crit) 
^2 

'-'O(crit) 


0.153 0.135 0.120 0.114 0.105 0.105 0.123 
0.150 0.141 0.132 0.129 0.123 0.123 0.135 









TABLE V: The critical relative shear square Eq^^^,|^j and 
the corresponding absolute shear square crQ(j.^i^j for Type-B 
bounce universe. If the initial shear exceeds this critical value, 
the universe evolves into singularity. 
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